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Abstract 

Diffusion coefficients are obtained from linear response functions and from the quan- 
tal fluctuation dissipation theorem. They are compared with the results of both the 
theory of hydrodynamic fluctuations by Landau and Lifschitz as well as the Boltzmann- 
Langevin theory. Sum rules related to conservation laws for total particle number, 
momentum and energy are demonstrated to hold true for fluctuations and diffusion 
coefficients in the quantum case. 

1 Introduction 



For a theoretical description of phenomena like nuclear multifragmentation mean field theories 
have to be extended to include the dynamics of fluctuations. In the case of the semi-classical 
Boltzmann-Uhling-Uhlenbeck (BUU) equation a first step was made by Bauer et al. [Q. Their 
numerical realization of the collision term corresponds to adding to the BUU collision term a 
further contribution. The average of the latter should vanish since BUU is meant to describe 
the effect of the collisions on average. In this sense the additional term acts like a fluctuating 
force. However, in jl] detailed properties of the former as its second moment (or correla- 
tion function) were not investigated. Moreover, in a later work the fluctuations that are 
obtained within the method of Q were judged unsatisfactory. 
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Possibilities to extend theories for average dynamics were proposed before by Landau and 
Lifschitz || for hydrodynamics, by Abrikosov and Khalatnikov jp for Landau Fermi Liquid 
theory and by Bixon and Zwanzig || for the Boltzmann equation. These authors augment 
the equations for the averages with a fluctuating term and derive explicitly its correlation 
function, albeit in different ways. 

The task to describe heavy ion collisions, especially multifragmentation, initiated the 
more recent approaches by Ayik and Gregoire ||, Randrup and Remaud ]7J and Hofmann, 
the author and Tsekhmistrenko || . In |]7| and || the equation of motion for the fluctuations 
(the second cumulants) of the distribution function (the Wigner transform of the one particle 
density matrix) is addressed, as first step. The inhomogeneity of this equation, called diffusion 
coefficient in analogy to Fokker Planck equations, can be related to the correlation function 
of the fluctuating force. The latter is derived in || by investigating the hierarchy of equations 
for the n-particle densities. It is found to coincide with the one obtained from the expression 
of J7| for the diffusion coefficient - for the case that the equation for the average distribution 
function is given by the BUU equation with a Markovian collision term. 

In the approach of [Q the diffusion coefficient is determined to guarantee that stationary 
solutions of the equations of motion for the dynamical fluctuations are, for a stable system, 
equal to the equilibrium values. The latter can be determined by means of the fluctuation 
dissipation theorem || to include quantum effects beyond Pauli blocking. In a recent work 
[TCfl these quantum effects were investigated in the density range between one tenth of up 



to two times the saturation density. They were found to be sizable for temperatures below 
15 MeV (This temperature range is of special interest as it is the one in which, in a certain 
density range below saturation density, homogeneous nuclear matter can become mechanically 
unstable against isothermal density variations, in other words, the temperature range of the 
isothermal spinodal region of nuclear matter.). 

In the approaches of || and |[7| - both now known under the name Boltzmann-Langevin 
theory - correlation function of the fluctuating force and diffusion coefficient, respectively, are 
determined as functionals of the distribution function. The latter is allowed to be the one of 
a non equilibrium state of the system. However, the numerical evaluation of these functionals 
turns out to be too time consuming as to presently allow for an application to a realistic (six 
dimensional) system. Therefore, recently JXTJ the expression for the diffusion coefficient was 
proposed to be evaluated at the equilibrium distribution, similar to the approach of ||. 

It is the main aim of the present paper to work out the relation of the approach of 
H| based on the (quantal) fluctuation dissipation theorem to those of Landau-Lifschitz and 
Bixon-Zwanzig as well as to the Boltzmann-Langevin theory. The connection to the latter 
will be established for an equilibrium state in the case that the classical form of the fluctuation 
dissipation theorem is applied. 
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In Boltzmann-Langevin theory sum rules reflecting conservation laws for total particle 
number, momentum and energy are fulfilled by the expressions obtained for the correlation 
function of the fluctuating force || and for the diffusion coefficients |T2|]. These sum rules 
will be demonstrated to hold true also for fluctuations and diffusion coefficient found in the 
approach of - in the classical limit as well as for the quantum corrections. A proper 
interpretation will be given for previous results || which seem to violate these sum rules. 



2 Review of general formalism 



In this section the basic ideas of the approach of || [Uj will be reviewed using a more gen- 
erally applicable formulation. It will be explained in which way equations of motion for the 
fluctuations are obtained starting from the equations for the first moments or averages. The 
former turn out to be inhomogeneous. The knowledge of this inhomogeneity, the diffusion 
coefficient, can serve as the first step in the determination of the correlation function of the 
fluctuating force. Moreover, the equations obtained can be used to study the behaviour of 
the dynamical fluctuations, at least in the range of small amplitudes |13| . 



The starting point are the equations of motion for averages A v of a set of dynamical 
variables A u , 

Mt) = m t , 

the system is to be described with. These equations are assumed to have the form 

d 

—A u + V VV ,A V , + c v = . (2) 

Here and in the following summation over repeated indices is understood to be performed. 
The linearity in A v may have been obtained by linearization. The c u are independent of A v 
and symbolize terms containing either external fields or zero order contributions remaining 
from the linearization procedure or both. 

To have a specific example in mind let us look at the case where the system is described by 
the distribution function (which is the Wigner transform of the one particle density matrix) 
np(r,t). In the semi classical limit the equation of motion of the latter is the BUU equation 
which is nonlinear in np(r,t). One, therefore, has to perform a linearization to obtain an 
equation that can be written in the form of (§). The averages A v would then be the deviation 
Snp(r,t) of the distribution function from the reference distribution chosen, e.g. the one of 
global equilibrium, 

A v {t) 6n p (r, t) = n p (r, t) - n°(e° p ) . (3) 

The index v stands for position r and momentum p and eventually for a spin-isospin index, 
the latter being suppressed throughout the paper. So far, the reference distribution might be 
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the one of a local equilibrium or some other state as well. One simply would obtain expressions 
for T> UfJ/ and c u different from those found when linearizing around global equilibrium. 
The fluctuations of the dynamical variables are defined as: 

a u ,(t) = (±[A u ,Al] + ^-A u (t)A;(t) (4) 

From this form the property 

°Vm = °\iv ( 5 ) 

follows even for non hermitian A u , as one has to deal with considering the spatial Fourier 
transform of 5np(r, t) . 

Since the equations of motion (fj) for the averages A v have been assumed linear, the 
corresponding ones for the A u differ from the former only by a term, called the fluctuating 
force, whose average vanishes. Thus the equations for the fluctuations are found to read 

9 r \ * 

— a VfJti + V vvl o v i^ + (V^a^ u ) = 2 d vft , (6) 

where the diffusion coefficient d Ufl includes all terms arising from the fluctuating force. In 
general, d Ufl depends on the momentary state of the system: looking at the hierarchy of 
n-particle densities mentioned before it becomes clear that d Ufl stands for terms containing 
the correlated part of the three particle density. To calculate these terms exactly would be 
equivalent to solving the full hierarchy which is, of course, not possible. One, therefore, has 
to use approximations. The one made in the present approach is the following. For a stable 
system d vli is chosen in such a way that the stationary solution of (|6|) is equal to the equilibrium 
fluctuations; for an unstable system d Vil is obtained [IIJ as the analytic continuation of the 



expressions found for the stable system. This is achieved by setting 

2d v , = V uv ,af,^(p^a%y (7) 

where, for the stable system, a^ 1 stands for the equilibrium fluctuations, while it represents 
the analytical continuation of the latter in the unstable case. 

The diffusion coefficient obtained in this way depends on both the equilibrium properties 
- via - and the dynamics of the system - via T> vtl . Let us look again at the example 
previously discussed where the linearized BUU-equation is the one from which V UfJi is to be 
obtained. In this case T> Vil containes information not only about the drift terms (involving 
the equilibrium mean field) but also about the collision term and thus about the differential 
scattering cross section. Therefore, the diffusion coefficient depends on both the mean field 
and the residual interaction causing the two particle collisions. 

Clearly, this choice for the diffusion coefficient is an approximation to the exact term 
appearing as inhomogeneity in (|6|). However, in the case that the system approaches equi- 
librium, the value of the latter approaches the one of the former, which thus can serve as an 
approximation valid for systems not too far from equilibrium. 
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Since T> u/1 is determined by the form (0) of the equation of motion for A u , what is left in 
order to calculate the diffusion coefficient is the determination of the equilibrium fluctuations 
(and their analytical continuation to the unstable case). For this purpose the relation of 
the latter to dissipative processes in the average dynamics as expressed in the fluctuation 
dissipation theorem is exploited. The latter provides the connection between the equilibrium 
fluctuations and the corresponding response function which can be obtained from the equation 
of motion for the averages A v as follows. 

Introducing an external field f/ ext in such a way that the change in the Hamiltonian due 
to U ext is given by 

5H = AlU^(t), (8) 

one defines the response function as the negative functional derivative of the Fourier transform 
A v \ (u) of the average A„(t) with respect to the external field: 



(9) 

c/ cxt =o 



8(A V ) (u>) 



This derivative can be calculated from the equation of motion of the A v . The dissipative 
processes are related to the quantity 

xlM = \^M-x»Mf) (io) 

which, therefore, is usually called, for real arguments u, the dissipative part of the response 
function (^). 

According to the fluctuation dissipation theorem || the a*^ are given by the sum of two 
contributions: 

< = < + A^, (11) 

to be explained in the following. 

The first contribution on the right hand side shall be given by 



where T stands for the temperature. The integral in (12) is to be taken along a contour C. 
For a stable system C lies on the real axis. In this case the stationary solution a st ( |TT] ) can 
be interpreted as the equilibrium values of the fluctuations and (|TTD with fll2]) is one version 
of the fluctuation dissipation theorem. For an unstable system, C has to be deformed in 
such a way that it lies above all poles of the response function x(o>) an d below all poles of 
the function xi 00 *) i n the complex frequency plane, crossing the line Re u = in the interval 
-2tt T/h < Imcu < 2-nT/h. 
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Let us turn now to the second contribution to cr st , i.e. to Acr^. Its physical origin lies 
in the possibility that, due to conserved quantities, an initially equilibrated system does not 
relax to the new thermal equilibrium (grand-, micro- or canonic) after an adiabatic switch on 
of an external perturbation. 

To shed some more light on this let us look at the case where the equilibrium distribution 
is assumed to be the one for the canonical ensemble. Then Aa UfJi is given jyfl by the positive 
semi definite matrix[] 

Aa Uf , = T( X ^-X^ = 0)) (13) 

where the isothermal susceptibility x T can be calculated from the behaviour of the equilibrium 
values of A, as functions of the external field: 



T _ ( 6<A V >^ 



(14) 

[/ext =0 



Here, the change in the equilibrium values of A u has to be taken under the subsidiary condi- 
tion of a constant temperature. The difference between isothermal suszeptibility and static 
response as appearing in flTj[ ) can then be expressed in the following way: using the isother- 
mal suszeptibility as scalar product in the space of dynamical variables one can determine 
an orthogonal basis {V c } in the subspace of those variables representing conserved quantities. 
One then finds fll5| : 

- x,> = o) = E ( 15 ) 

c Xv c v c 

where the sum on the right hand side extends over all elements V c of the basis. This formula 
clearly states the relation of Act to the existence of conserved quantities. 

These considerations are valid for a system in contact with a heat bath preserving the 
temperature. However, if the system to be investigated has no heat bath, one might conclude 



that, for such a case, Aov M has to be modified as compared to (p~3|) , like it was done in |L0 . 

In the present paper the following point of view is taken: Quite generally a contribution 
Act to CT st must be expected to exist; there will be no calculation of the former; for the 
comparison with the other approaches only <j st will be used. It will be found that the former 
are consistent with Act = 0. Actually, we would like to use this result as justification to 
identify CT st (|T2|) with a st (|H]). Nevertheless, a st will be distinguished from <r st by keeping the 
tilde whenever the expression on the right hand side of fll2|) is to be referred to. 

At this stage it is possible to briefly compare with the approach of Bixon and Zwanzig 
[[J. (For more details see the appendix A.) These authors address a stable system with the 
Boltzmann equation as the one governing the average dynamics of the distribution function. 
For comparison of the results of the two methods one, therefore, has to define T> Ufl and c„ by 

1 A matrix M is called positive semi definite, if J2 V u My^XyXp > f° r anv rea l x v 
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writing the linearized Boltzmann equation in the form (0) with (|3]). One finds the product of 
the expression on the right hand side of ([7|) with a 5— function in time to equal the expression of 
[[| for the correlation function of the fluctuating force. This is consistent with the assumption 
of a vanishing memory for the fluctuating force made in |Q. The value for the equilibrium 
fluctuations used in || is the same as the one obtained from flP2|), provided the classical limit 
of the latter is taken and the response function is calculated from the Boltzmann equation 
(Please notice, that in the form of the latter used in || no particle interaction is taken into 
account.). This can be concluded from the results in [|lO|j . 



3 Fluctuations for Hydrodynamics 

In this section the method described in the previous one for the determination of the diffusion 
coefficients is applied to hydrodynamics. Agreement with the method of Landau and Lifschitz 
[[| will be established. For the case of hydrodynamics the index v stands for a continuous 
variable and a discrete one. The former is given by the spatial position r or - in the corre- 
sponding Fourier space - the wave vector k. The discrete one is meant to tell whether A v is 
the deviation bp of the mass density from the equilibrium value po, the deviation 5q = T p 6s 
of the density of heat or the momentum density g. (Here s is the entropy per unit mass and 
To is the equilibrium value of the temperature.) 

The linearized hydrodynamic equations for the dynamical variables, i.e. those correspond- 
ing to A u , read (see e.g. §132 of [0]): 

^8p(k,t)+ikg(k,t) = (16) 

d 

— gi (k, t) + i kj Uy (k, t) = i kj Sij(k, t) (17) 
d 

— 5q (k, t) + i kj T (k, t) = -i k G(k, t) (18) 
Here beside the usual stress tensor 

Ikj(k, t) = 6y P(k, t) - t(D t (% g, + h gj ) + (A - 2 A) Sy k g) (19) 

and heat current 

j T (k,t) = -ikK5T(k,t) , (20) 

where P is the pressure, Di and D t are the longitudinal and transversal diffusion constants, 
respectively, and ST is the deviation of the temperature from T , the fluctuating forces Sij 
and G appear. In the form suitable for later purposes their correlation functions are given 
iby: 

\[s i3 {k,u) , s mn {-k',t)] + \ = (2rrf5(k - k') 2E{u) e iujt p D ijmn (21) 
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Gi(k,u) , Gj(-k',t) , > = (27r) 3 (5(fc - k') 2E(u)e iuJt KT 5 lj (22) 



1 

2 

where 2£(u;) = hu coth(huj / 2T) and Ajmn = D t (5 im 5 jn + 5 in 5 jm ) + (A - 2D t )SijS mn . The 
mixed correlation vanishes: < G Sij >= 0. In || the correlation functions are obtained using 
their relation (e.g. §121 of |1(J) to coefficients appearing in the expression for the rate of 
change of the entropy (see §119 of jl6| ). 

The appearance of the factor E(u) implies that in the quantum mechanical case, i.e. for U ^ 
0, the correlation function exhibits, in general, a finite correlation time. Furthermore, E(u) 
reminds one that diffusion constants Di, D t as well as thermal conductivity k are functions of 
frequency, the even parts of which decrease for large frequencies asymptotically, at least, as 
1/uj (for the right hand sides of (|21|) and ( p2]) to be proper Fourier transforms). This feature 
ensures convergence of frequency integrals which will have to be dealt with in the following. 

The structure of the equations for the averages follow by putting the right hand side of 
(|16l)-(|rS|) equal to zero. The equations of motion for the fluctuations of dp, g and 5q can 
be obtained in different ways. One of them is to directly use the equations (ffBI)-(ffq) for the 



dynamical variables exploiting the relations (pTJ) , (22) for the fluctuating forces. This one will 
be applied later. An other possibility is the one described in section 2. For the latter one 
needs to know the corresponding response functions, as the density density response. 
If external fields are introduced in the equations of motion according to the rules 



p + poU; 

g^ g + p U c (23) 
5T^ 5T + Tq 



then the interaction energy takes the form ||17|| : 



H ext = Jdr (Sp U e p + g U c + 5q Ufj (24) 

This is exactly of form (|8|). Therefore, the equations for the response functions can be obtained 
taking the functional derivatives of (|16|)-(|T8|) with respect to the external fields introduced 
according to d23|) . For example, the set of equations containing the density density response 
X PP (k, k',u) = (2n) 3 5(k — k')xpp(k,u) is found by calculating the derivative with respect to 
Up. The solution reads: 

-k 2 

Xpp(k,uj)=p 1— (25) 

PP LU 2 + tLuD l k 2 -C 2 k 2 + tdk^C 2 ^^ 

if the following abbreviations are used for derivatives of thermodynamic functions: 

k T (dS\ 2 (dP\ (dT\ , , 
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The other response functions are found to be related to xpp i n the following way: 

Xg iP (k, u) = x Pgl (k, u) = ki — x PP (k, w) (27) 

X qp (k, w) = X Pq (k, w) = - - . dk<2 Xp P (k, u) (28) 

X qgl (k, u) = X Si q(k, uj) = ki^ Xq P {k, u) (29) 
—inak 2 inT k 2 

X ^ u) = ^TdT 2 x ^ ^ + ^TidV (30) 

X Sigj (k,u) = j^xi(k,uj) + ^ - Xt(k,uj) (31) 

U) 2 

Xi(k,uj) = — Xpp(k,u) + po (32) 

(U \ Wtk2 /QQ\ 

From these functions the equilibrium values of the corresponding fluctuations can be calcu- 
lated - at least the contributions a st - combinations of which form the diffusion coeffi- 
cients. 

According to ©, with the operator V defined by (||) and (|ilf )-(|ilj|), the equations of 
motion for the fluctuations read: 

d 

— a pp (k, k) + i kj a gjP - i k j a pg] = (34) 

d 

—a pgi (k,k)+i kj a gjgi - i fcj a pUtj = 2 d pgi (35) 



dt 
c 

dt 



d , 2 

—<r pq (k, k) + i kj a g]q + k k a pT = 2 d pq (36) 



d 

—a gigj (k, k) + i h a Uagj - i k[ a giUjl = 2 d Sig . (37) 
d 

■7^<Tg iq (k, k') + % kj a Uijq + k k' 2 a giT = 2 d gtq (38) 

9 2 

—a qq (k,k') + Kk 2 a Tq + Kk' a qT = 2d qq (39) 

To simplify notation the arguments k, k' are written only in those terms containing a deriva- 
tive with respect to time. Please notice that in the present case the relation @ translates 
to <jAB(k,k') = aBA(—k',—k), where A, B = p,g,q. To shorten the notation further the 
abbreviations 

h crn a A{k, k') = ki cr PA (k, k') - iD imjn k m k n a SjA {k, k') , (40) 
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2 , PO 

opa = c a p A + — aa qA 

-*0 



and 



a at = a&Ap + 



1 



Poc v 



O'Aq 



(41) 
(42) 



2 d pgi (k,k) — ikj cr g . gi i kj — i kj a^ gi i k i afp 



are introduced for certain combinations of fluctuations. 

The diffusion coefficients follow either from the general formula (|?]) together with @, 
"l-(ffED or from (M]-(M): 



(43) 

(44) 
(45) 
(46) 
(47) 



2 ( ; k ^ — X (Tg . ^ — I - /x ^ pT — ^ ^ ^ pT 

2 °^gigj {k,k) = ik\ On^gj — & ^giix,-; = Di m i n k m k n a gig . 

2 d giq (k, k) = ikj <J^.. q + K,k a s giT = iki a s p q 



B, ■ k' k' rr st 



m-"n" gig ; 



2 d qq (k, k ) 



«: fc 2 + k k' 2 



> Tq -rrvnu u qT 

where a s ^ B = o- s ^ B {k, k'). To obtain these results use was made of the behaviour of p,q,g 
under time reversal. 

As mentioned earlier in this section, knowing the equations (|T6"l) -(|l8l) for the dynamical 
variables one can determine the equations for the fluctuations and the diffusion coefficients 
in a second way. Realizing that the time dependence of o VVL (t) can be interpreted as average 
of time dependent dynamical variables A u (t) one can calculate the time derivatives of the 
fluctuations by inserting the expressions for time derivatives of dynamical variables. The 
resulting equations for the fluctuations are identical to (j34])-(|3~9l). For the diffusion coefficients 
one finds the following explicit formulas: 



2d h pgi (k,k') 



Sp(k,t) , Sij(-k',t) 



2 d L (k, k') = i k'- ( \dp(k, t) , Gji-k', t) 



2d L (k,k')=tk l - 



2d L giq (k,k') 



2 

3 2 



s u (k,t) , gj(-k',t) 
Sij(k,t) , q{-k',t) 



2 d L qq (k, k') = -i k 3 - ( Gj(k, t) , q(-k', t) 





(48) 




(49) 


-k',t) 


+ ) ( 5 °) 


-k',t) 


+ > (51) 


-k',t) 


L) ( 52 ) 



To evaluate them one has to express 5p, g, Sq as functions of the fluctuating forces with the 
help of (pT6|) (pT8|) . Using the response functions (|25|), (|2?|)-(|3"3"D one can summarize the result 
by 

A(k, to) = XA P {k, uS) — ^ Sij(k, u) + XA q {k, Gi(k,u), (53) 



Po k 2 



kT k 2 
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where A — p, kg, g, and 



00 



Xt{k,uo) 



a; 



(54) 



d L AB (k,k') = {2uf5{k-k') d AB {k) 

where A,B — p, g, q. For the d(k) one ends up with the following expressions: 

duo 



d Pgl ( k ) 

dpq(k) 

d qq (k) 
di(k) 
d t (k) 



Using (0), (|54]) as well as the correlation functions ( plj) and ( ^2|) one finds that the diffusion 
coefficients (HS|)-(B2p are diagonal in k, k': 



(55) 

(56) 
(57) 
(58) 
(59) 
(60) 
(61) 



o E(lo) (-iDi) k { Xppik, oo) 

duo . . 

— E(w)Xpg(fc,^) 

cL; uo + iDik 2 

2k ^ w Xgl9 ^ ' U > 

duo 

— 2E(oo)(-iD l ) ki x glP (k, 



2E(u)x qq {k, u) 



duo 
2n~ 



2E(Lu) X t(k,Lo) 



With the help of the expressions (|25|) , ( 
results for the diffusion coefficients (^3| 
is inserted, i.e. 

d AB (k, k') = (2n) 3 5(k - k') d AB (k) = d L AB (k, k' 



j) for the response functions one finds the same 
provided that for cr st only the contribution a st 



(62) 



with A,B = p, g, q, where the tilde above d reminds of the restriction on a st . 

We thus see, that the theory of fluctuations of Landau and Lifschitz leads to the same 
expressions for the diffusion coefficients as our method, provided the additional contribution 
Aa to <7 st ( ITT]) either vanishes exactly or is neglected. 



4 Comparison with Boltzmann-Langevin theory 

In this section connection will be established between the approach explained in section 2 and 
the Boltzmann-Langevin theory by comparing the expressions for the diffusion coefficients. 
To calculate the latter with the method of section 2, in general one has to know the behaviour 



of the response function as function of frequency. However, investigating (]T2|) one finds that 
less information is needed to exploit the fluctuation dissipation theorem in the classical limit: 
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one only has to determine the static response, i.e. the value of the response function at u — 0. 
From (|J) and @ the equation for the latter follows as: 

5c 

?W> = 0) = ^ xt( J =Q) (63) 

The combination on the left hand side of (|63"D appears also in the expression (]?]) for the 
diffusion coefficient since the contribution 5~ st ([12]) to the latter equals, in the classical limit, 
the product of static response and temperature. Therefore, it is not necessary to solve (|63|) 
for the static response - only the right hand side of (|63|) is needed. 

Let us now look at the cases of BUU and Landau equation as the ones for the averages, 
i.e. the ones defining c v via (Q). Linearizing the equation for the distribution function in the 
deviation from the homogeneous equilibrium distribution one finds 

— 6n p (r, t) + (v p V r - Jp,r*) \Sn p (r, t) - g^Se p (r, t) \ = . (64) 
The linear operator J is defined by the collision term: 

I[n p ] = Jp, r * (n p (r, t) - n°(e p fj (65) 

The asterix reminds that, in general, Jp r includes integrations over momentum and position. 
(In the non-Markov case, there will be an additional integration over time. An index t at J, 
however, is omitted for the moment.) 

Please notice that the energy appearing in the argument of the equilibrium distribution 
n° is taken to be the momentary one. The latter differs from e° which appears in @ and 
which represents the energy in equilibrium without external field, by the amount 

Sep = J j^fpp'Snp, + U p xt . (66) 

The kernel fppi is given by fpp = dll/dp for the case of the Boltzmann equation with a 
momentum independent mean field U[p\. In Landau Fermi Liquid theory fpp> is called quasi 
particle interaction. 

To calculate the diffusion coefficient one has to determine c v as functional of the external 
field. Comparing (|M|) with the general form (^|) of the equation for the averages one obtains 

v p V r - Jp, r *)[-^- p U p( r > *)) • ( 67 ) 
The right hand side of the equation (|63|) for the static response corresponds therefore to 

m J^- W) « (68) 
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Denoting by d the contribution from a st ( |T2"D to the diffusion coefficient (|7|), one finds from 



2d(p,r,p',r')\ ci = -T (j p , r * (-^h 3 6(p - p')8(r - r')j + (r,p 




(69) 



1 f dp' 

Jp,r * h(p, r, t) = -— h(p, r,t) + J -j^Ipp' Kp » r > *) 



(70) 



since the terms with the spatial derivative in (|6q) cancel because of (V r + V r ')5(r — r') = 0. 

For a local, Markov collision term the linear operator J Q55| ) can be written as the sum of 
a diagonal and a non-diagonal contribution 

dp' 

ayp, r, i) -f / 

where /i(p, r, t) is some function of p, r, £. The action of the diagonal part of J consists of a 
multiplication with a function of p only which, in ([fT|) , is written as 1/rp to indicate, that 
this factor has the dimension of an inverse time. It is uniquely determined by J under the 
subsidiary condition that the kernel lpp> does not contain terms proportional to 8(p — p'). 
For the collision term fl70|) the diffusion coefficient follows as: 



d(p 1 ,r 1 ,p 2 ,r 2 ) = S(rt - r 2 



Pi 



P2) - 2 ( J PiP 2 n S 2 < 



+ ip 2 pAiKi 



(71) 

using the definition Hp = 1 — n p . 

This result can be compared with the one obtained in the Boltzmann-Langevin theory for 
the (Markovian) BUU collision term. The comparison will be performed in two steps: In a 



first step results obtained in |11J for the low temperature limit are used. In a second step the 
restriction to small temperatures will be removed. 



In [|TTJ the Boltzmann-Langevin expressions for collision term and diffusion coefficient were 
evaluated in the low temperature limit and for small deviations from thermal equilibrium. The 
result is written in formula (21) of [|TT|. To perform the comparison one first has to determine 
Tp and lpp> as defined by flTOl) . Comparing |70| with the expression for the collision term in 
p| (there called drift coefficient) one finds that, in the approximation of fil l, 1/rp is given 



by the momentum independent quantity W of [|TTJ and the kernel I Pl p 2 is given by A°/i^i 2 
of ||11|| . With the help of these relations the expression (|7lD turns into the lower one in formula 



(21) of [11], i.e. into the Boltzmann-Langevin result for the diffusion coefficient in equilibrium 
and in the low temperature limit. 

However, the restriction to low temperature can be released. To this end one starts from 
the general Boltzmann-Langevin expressions for drift- and diffusion coefficient Hl2 , §]. One 
finds that the contribution to the collision term which is linear in 5np can be written as (see 
appendix B) 



I\n 



pjllin 



a 2 [n°](p] 
2n° p n° p 



5n p - 



dp 2 Ocojnj (P,P 2 ] 
h 3 2n° 2 n£ 2 



5n 



P2 



(72) 
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where (in the notation of [|12|]) a[n](p,r,p' ,r') = 5(r — r')(a 2 [n](p) h 3 5(p—p')+a cov [n}(p,p')^ 
appears as inhomogeneity in the equation for the fluctuations, i.e. is the equivalent to 2d(p,p') 
of the present paper. Comparing (ff2|) with (|70 ) one finds the identification: 

"») 1 f H</t P 

Boltzmann— present 

T ~ V (?3) 

Langevm work 

Inserting (|73|) in the expression (ffll) results in 

2rf(p,r,p , ,r / )|ci = a[n°](p,r,i/,r') (74) 

one only has to use the symmetry property a cov [n°](p,p') = a cov [n°](p' , p) . Thus, the equi- 
librium diffusion coefficient of the Boltzmann-Langevin theory turns out identical to the 
contribution d to the one of the present approach provided the latter is calculated with the 
classical form of the fluctuation dissipation theorem, i.e. is identical to ([H|). 

As mentioned before, in the Boltzmann-Langevin theory the correlation function for the 
fluctuating force, C(p,r,t,p' ,r f ,t r ), when calculated in equilibrium and for a Markovian 
collision term, is given || |7|] by 

C(p, r, t, p', r', = 2 d(p, r, p' , r') S(t - t') , (75) 

cl 

where the relation between the Boltzmann-Langevin diffusion coefficient and (|7l| ) was utilized. 
Using the method of Landau and Lifschitz, applied in the previous section, to calculate the 
correlation function of the fluctuating force for BUU or Landau equation (for the latter see 
e.g. one obtains a result that agrees with (|75f ), again for the classical fluctuation dissipation 
theorem and for a Markovian collision term. 



In a recent work |L8| the Boltzmann-Langevin theory is extended to a non-Markovian 
collision term. In the framework of the theory of Landau and Lifschitz this generalization can 
be done without difficulty, since until the final step in the derivation one does not need to 
specify the form of the collision term. Indeed, one finds for C(p,r,t,p' ,r' ,t') an expression 
similar to the right hand side of fl69l) . However, in the former an 5-function S(t — t') appears 
besides those in the positions and in the momenta, and the operator J fl6"5D is now non local 
in time: r / r 

J p>r ,t * h(p, r,t)=J P ^ 1 l(p, r, t; p', r', t') h(p', r', f) (76) 

Following || the quantum correlation function of the fluctuating force is obtained multiplying 
the Fourier transform (in time) of the classical one by E(u)/T. Thus one finds the former to 
read: 



C(p,r,t,p',r',t') = J !te-**E(u,) J dre™ Jp, r , T * {^h 3 5(p - p')S(r - r')S(r - t>) 



+ 



+ (p,r,t^p',r',t') (77) 
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In the expression for C the use of the kernel I ( |76|) is omitted for the sake of an easier 
comparison with (|69|) in the Markov limit. 

It is straightforward to show that (|77D leads to the result found in |L8| for the correlation 
function of the Boltzmann-Langevin theory in equilibrium, provided the expression for the 
non Markovian collision term of |18) is used to define J according to (|65|). 



5 Sum rules — the case of zero wave vector 

In finite, closed systems the total number of particles, the total momentum and the total 
energy are conserved quantities. Starting with an ensemble of systems all having the same 
particle number, momentum and energy, the dynamical variables corresponding to the latter 
should not exhibit any fluctuations at a later time. This feature constrains the possible forms 
for the equations of motion of the fluctuations. Especially, it implies sum rules for fluctuations 
a and diffusion coefficient d. Since the deviations of total particle number, momentum or 
energy from the corresponding equilibrium values are given by the values of the phase space 
integral of the deviation of the distribution function weighted with 1, p or Ep, respectively, 
the sum rules read: 

drdp 



J^f-O p b(p,r,p,r')=0 (7c 



where Op = l,p,ep and b = a,d. In the following it will be demonstrated that the results 
for a and d obtained with the method of section 2 together with the assumption Act = are 
consistent with these sum rules. 

At first, the classical limit shall be addressed. Here, the diffusion coefficient obtained in 
section 4 was found to agree with the result of the Boltzmann-Langevin theory in equilibrium. 
The former is, therefore, known |], to fulfil the sum rules (|78|). 



For investigation of the fluctuations one might think of solving (p3|) for the static response. 
However, there is some subtlety involved: response functions may be non analytic at uj = = 
k and one has to choose carefully the correct sequence in performing limits. This fact is 



relevant for the present case since in the sum rules (|78|) there occurs an integral over space. 
The value of the latter is thus equal to the one of the Fourier transform of the integrand taken 
at k = 0. Moreover, in the contribution cr st to the integrand there appears, in the classical 
limit, the response function at uj = 0. The order in which the two limits - k, u — > - have to 
be taken can be found as follows. As said before, the sum rules arise due to the conservation 
laws for particle number, momentum and energy. The deviations of the latter from their 
equilibrium values are the averages of dynamical variables given by 

M Op (k = 0)= J^-Op5h p (r) , (79) 
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where Op = l,p,ep respectively and < 5hp(r) >= 5np(r). Due to the conservation laws 
M Qp {k = 0) does not exhibit any fluctuations: 

o~M 0p (k=o),A = , (80) 
where A is some dynamical variable. Inserting here (79) one recovers (f78|) for b = a when A 



is taken as np. The fluctuations are related to the static response whose microscopic form is 



given by 1 17 



XM Op (k= )A" = o) = fdt % - ({M 0p (k = 0) , A\-t)\) S(t) 



n 

The expression on the right hand side can be rewritten to give 

drdp lkl 



XA/ Op (fc=0),A(^ = 0) 




9(f) . (82) 



If the limit k — > can be commuted both with the averaging procedure and with the time 
integration, then one finds: 

XM Op (k=0),p'(u = 0) = lim ^lim J ^ O p xpp>(k, u^j (83) 

Therefore, the sum rules for the fluctuations (b = a st ) read in the classical limit: 

Acr MOp (fc= o),p' + T lim Him J ^ Op\pp (k,u)J =0 (84) 



Thus, for calculating the left hand side of (|78|) one has first to perform the momentum integral 
of the response function at finite k, uj weighted with Op. Then one has to take the limit k —>■ 
(corresponding to the spatial integral) and, finally, u —>■ 0. Due to this sequence of limits the 
investigation about the fluctuations fulfilling the sum rules turns out to be not easier in the 
classical limit than in the general case. 

There is a further, indirect argument that the sum rules are related to the order of limits 
stated, i.e. first k — > and then uj — » 0. It is well known ||19|| , that taking first the limit for 



uj and then for k leads to thermodynamic derivatives, e.g. to the compressibility in the case 
of the density density response function. However, a finite, closed system can have a finite 
compressibility, too. And indeed, solving ( |63D in the case of Landau Fermi Liquid theory one 
obtains M the known, finite value for the compressibility pCfl . 

In the general (i.e. not classical) case one may proceed as follows in order to prove the 
sum rules to hold true for fluctuations and diffusion coefficient. Let again Op be either 1, p 
or Ep. Multiplying the equation of motion for the distribution function rip(r,t) by Op and 
integrating with respect to r,p one finds [^(J 

d -^Op§- t 5n p (r) = 0. (85) 
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To obtain this result one only has to integrate by parts the terms containing the derivatives 
with respect to r and p and to use the feature of the collision term that number of particles, 
momentum and energy are preserved. 

Taking the functional derivative of (|S5D with respect to the external field results in 



§if^r0pxpp'(r,ty,t') = 0. (86) 



Since the response function is proportional to the Heavyside-0-function, Xpp'( r ^: r '^') ~ 
0(i — t') [[T7|], reflecting causality, the solution of ([86]) reads 

^O p xpp>(r,t,r',t') = 0. (87) 



Due to ([T0| ) the dissipative part of x has the same property. Therefore, it follows from (|T^) 
that 

f^r°P S PP<r> rr ) = 1% O p a^ pl (k = 0) = 0. (88) 

We thus see that the contribution <r st , when being calculated from the fluctuation dissipation 
theorem, fulfills the same sum rules as the fluctuations themselves. An approximation like 
Act = 0, in other words approximating cr st by CT st , is, therefore, consistent with the sum rules 
( |78D for the fluctuations. 

Let us now turn to the diffusion coefficient. According to (0) the latter is given by 

vpp>(r,r') - -Q^fpp"(r,r")appr(r",r')j + 

+ (p,r «-> p',r'Y (89) 

where integration over double primed quantities is to be understood. The last term on the 
right hand side stands for the complex conjugate of the first term but with the unprimed 
variables exchanged with the single primed ones. Multiplying d with Op and integrating the 
product with respect to position and momentum one finds that the contribution from the first 



term on the right hand side of (89) vanishes. This is due to the same reasons that lead to 
(|S1)D, i.e. directly due to the structure of the equation of motion for the distribution function. 
The other contribution vanishes due to ( p8|) with (|5|), i.e. due to features of the fluctuations. 
Thus the sum rules ([78]) for b = d 

J^rO p d(p,r,p,r') = J^O p d(p,p',k = 0) = 0. (90) 

are proven. 

The diffusion coefficient forms the inhomogeneity in the equation for the (non equilibrium) 
fluctuations. Due to ( p0| ) the equation of motion for the product of the fluctuations with Op 
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integrated over r, p is homogeneous. Therefore, the dynamic fluctuations obey the same sum 
rules as a st , i.e. (|8~8j) , for all times provided they did in the initial state. 

It should be mentioned that these properties, especially the ones of the fluctuations, are 
restricted to k — 0. Indeed, it turns out [[13], that the limit k — > of o~pp>(k) differs from the 
value at k = 0. This means, that the fluctuations are not continuous functions at vanishing 
wave vector. 



6 Summary 

The problem to determine equations of motion of fluctuations and diffusion coefficients was 



addressed. A method || [L0| for its solution using knowledge only of the average dynamics 
was reviewed using a more general formulation. The basic ingredient of the method is the 
fluctuation dissipation theorem and a suitable extension of the latter to instabilities. This 
extension consists in the prescription to determine the diffusion coefficients as the analytic 
continuation of the expression valid for the stable system. 

The connection between the diffusion coefficient in the present approach and the corre- 
lation function of the fluctuating force in the one of Bixon and Zwanzig || was explained. 
The method was then applied to hydrodynamics. Equation of motion of the fluctuations and 
diffusion coefficients were demonstrated to be the same as those obtained using the theory of 
fluctuations by Landau and Lifschitz Subsequently, the case of BUU and Landau equation 
as equation for the average one (quasi) particle distribution was addressed. An expression for 
the diffusion coefficient in terms of the linearized, Markov collision operator was calculated 
using the classical form of the fluctuation dissipation theorem. This expression was found to 
coincide, for BUU, with the one obtained for the equilibrium distribution in the Boltzmann- 
Langevin theory jy]]. The connection with both the theory by Landau-Lifschitz and the 



Boltzmann-Langevin theory holds under the assumption that for the equilibrium fluctuations 
- and, consequently, their analytical continuation - only that contribution is taken into ac- 
count which arises from the frequency integral of the dissipative part of the response function 
(weighted with Ti coth(hu/2T)/ir). 

In the framework of the theory of Landau and Lifschitz the generalization to a non- 
Markovian collision term in Boltzmann and Landau equation was presented. The correlation 
function of the fluctuating force obtained for the case of the Boltzmann equation coincides 
with the equilibrium one of the extension [|ll| of the Boltzmann-Langevin theory provided 



the same form is used for the collision term. 

To establish sum rules arising from conservation of particle number and momentum the 
behaviour of the response functions for vanishing wave vector and frequency was taken into 
account properly. Clarifying the sequence in which the limits k —>■ and uj —>■ have to be 
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calculated, results of previous work || could be interpreted correctly. 
Acknowledgments 

The author wants to thank H. Hofmann for many discussions during the past years, for 
helpful criticism and suggestions as well as for carefully reading the manuscript. This work 
was supported in part by the Deutsche Forschungsgemeinschaft as well as by the National 
Science Foundation under the Grant No. PY-9403666. 



19 



Appendix: 



A Detailed comparison with Bixon and Zwanzig 

In this appendix the comparison of the present approach with the one of Bixon and Zwanzig 
[H summarized at the end of section 2 will be worked out in more detail. 

Bixon and Zwanzig start with the classical Boltzmann equation without mean and/or 
external field as the equation for the average one particle distribution function. Augmenting 
the linearized Boltzmann equation by a fluctuating force, fs F(t) in their notation, the authors 
write the result in the form 

^ + L<p = F(t) (A.l) 

where <j)(t) is the relative deviation of the distribution from the Boltzmann distribution. The 
correlation function of the fluctuating force is written as 

{F{k u v x ; ti)F(fc 2> v 2 ; t 2 )> = 2B{k u v 1; k 2 , v 2 )6(t 1 - t 2 ) . (A.2) 

For B the authors find the expression 

2B{k u vi, k 2 , v 2 ) = (La + L 2 ) (fafa) , (A.3) 

where the average on the right hand side has to be performed in equilibrium. For the calcu- 
lation of the correlation function (A.2) the authors use the expression 



(<Kfci, vO0(fe 2 , v 2 )> = [/ B (v 2 )]-^( Vl - v 2 )5(k 1 + k 2 ) (A.4) 

After this repetition the comparison can be performed easily. This the more, as the 
equation for the average distribution following from ( |A.1| ) by averaging of both sides, is 
written directly in the form (^): one only has to identify A v with cj)(k v ,v v ). Clearly, one has 

V^A„ <-> (L<f>)(k v ,y v ) . (A.5) 

From (0), ( |OD and (|AT5|) follows 

d vli = B(k u , v u , -kfj,, v^) , (A.6) 

where the minus sign in front of fc M is due to the hermitian conjugation of the second dynamical 
variable in the definition (||) of the fluctuations. Thus the statement at the end of section 2 
about the relation of the diffusion coefficient and the correlation function of the fluctuating 
force is proven. 



In the rest of this appendix it will be demonstrated how the expression ( |A.4j ) for the 



equilibrium fluctuations in a classical ideal gas can be obtained from the related response 



20 



function. First it should be noted, that in ( |A.4j ) there appears the fluctuations at finite wave- 
vectors ki,k 2 . Next one realizes that in the classical limit only the static response is needed 
to obtain a. From (|64j) one finds after Fourier transformation of the time and putting u = 0: 



5n p {r, u = 0) - — 5e p {r, u = 0) J = (A.7) 



The solution of this equation reads 



8n p (r, to = 0) - — Se p (r, to = 0) = c p (r) (A.i 



de p 

where (v r V r — J P) r*) c p( r ) — 0- This seems to be no progress at first site, but it is: J does 
not depend on the external field; consequently the same is valid for c p (r). Therefore, taking 
the functional derivative of ( |A.8| ) with respect to U ext (r',uj = 0) one finds the static response 
of non interacting particles after Fourier transforming the spatial dependencies: 

Xpp ik, k', u = 0) = -Q^hHtp - p')(27r) 3 5(k - k') (A.9) 



v 



Since the derivative of the Boltzmann distribution with respect to the energy equals the 
negative of the distribution divided by the temperature one ends up with the result for the 
equilibrium fluctuations: 

a(p,k,p',k')U id . gas = n°(p)h 3 5(p-p')(2n) 3 5(k-k') (A.10) 

where n°(p) is the equilibrium distribution function for the momentum with the property 
/ d 3 p/h 3 n°(p) = J d 3 vfB(v). Therefore, = (m/h) 3 n°(mv) and from ( |A.10| ) 

(Sf{yx, fci)5/(v 2 , fe 2 )) = f B (v 1 )S{v 1 - v 2 )(2n) 3 6(k 1 + fe 2 ) (A.ll) 

Dividing both sides by /b(vi)/b(v2) one finds ( |A.4j ), beside a factor (27r) 3 which is due to the 
different conventions of the Fourier transform. 
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B Collision term in BL-theory 



In this appendix the derivation of (|72|) will be presented. The notation used is the one of [0 . 

In a first step one observes that the linearized BUU collision term can be written in the 
form: 



I[n p ]\nn = f (t[dpi)^^(5{p-p 3 )-5{p-p 1 )^(Sn 1 nl+5n 2 n^ 



i=2 



(B.12) 

where as short hand notation a number i stands for p { and, as before, n p = 1 — n p . Using the 
quantities W x introduced in [|1^] one finds: 

I[np]W= -(w + [nl} + W-[n° p })5n p 

- J dp 2 [n° p (W ++ (p, 2) - W+-(p, 2)) + n° p (w -(p, 2) - W + (p, 2))] 5rip 2 , 

(B.13) 

where the index indicates, that the corresponding quantity is to be taken at n°. Making 
use of the properties of the Fermi function n° as well as of the energy conservation contained 
in the transition rate w one can rewrite the integrand in (|B.13|) as 



n°(W ++ (p, 2)-W+-(p, 2))+n°JW -(P, 2)-W + (p, 2)) 



n° p Zo(p,2) + n G p Z+(p,2) 
2n° 2 n° 2 



(B.14) 



Again the definitions of the Z ± are to be found in \Y1\. 

Remembering that the non diagonal part to a[n](p,p') is related to the via 

a cov [n°}(p,p') = n p Z (p,p') + n° p Z+(p,p') 

while the diagonal part is given by 

a 2 [n ](p) = 2ny p (w^n p ) + W-[n° p ]) 

one obtains fl72| ) from ( p,13|) . 



(B.15) 



^B.16) 
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